Understanding the Nucleation Mechanisms of Carbon Nanotubes in Catalytic 

Chemical Vapor Deposition. 
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The nucleation of carbon caps on small nickel clusters is studied using a tight binding model 
coupled to grand canonical Monte Carlo simulations. It takes place in a well defined carbon chemical 
potential range, when a critical concentration of surface carbon atoms is reached. The solubility of 
carbon in the outermost Ni layers, that depends on the initial, crystalline or disordered, state of the 
catalyst and on the thermodynamic conditions, is therefore a key quantity to control the nucleation. 
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Although multiwall carbon nanotubes are now pro- 
duced on an industrial scale [l| , a detailed understanding 
of the synthesis mechanisms of single wall carbon nan- 
otubes (SWNT) is still mandatory for more elaborate 
applications [2| . Diameter analysis of the tubes produced 
by Chemical Vapor Deposition (CVD) Q and a tendency 
towards chiral selectivity 0, 0] have been reported, but 
simple questions such as the role and physical state 
of the catalyst are still not answered. Transmission 
Electron Microscopy (TEM) observations performed 
after the synthesis lead to contradictory results showing 
evidences for either crystalline Q or liquid state 0] of 
the metal particle during the synthesis. Remarkable 
progress has been made in the in situ observation of the 
growth 0, 0, [h], 11 1, but the atomic resolution is not 



yet obtained under these conditions. 
Computer simulation is then a unique tool to gain an 
insight at this atomic scale that is very difficult to 
access experimentally. However, obtaining an accurate 
description of the interatomic interactions, as typically 
provided by first principles calculations, and addressing 
the size and time scales relevant to the experiments are 
extremely difficult challenges. Efficient catalyst particles 
for CVD are in the 1-10 nm diameter range and growth 
rates are in the nano- to micrometer per second range 
@; Such scales are largely beyond the capabilities 
of Molecular Dynamics (MD) techniques, meaning that 
only some elementary steps can be studied, or that the 
growth conditions imposed in the simulations are orders 
of magnitude too fast. 

Static ab initio calculations were used to study the 
carbon/catalyst interaction FJ, Il4 , 1 1 511 - Hofmann et 
al. [16| and Abild-Pedersen et al. [Fo| studied surface 
or subsurface C diffusion on Ni, and schematic growth 
model for CVD were proposed [3, [TtJ ■ First principles 
MD simulations were used to study the root growth on 
a Co catalyst [lU and the nucleation of a C cap on a 
Fe cluster was described by Raty et al. The latter 



explained the absence of diffusion inside the Fe particle 
by its nano size (55 atoms). However Ding et al. f2fij ]. 
using a semi empirical FeC interaction model, found that 
nucleating C islands on a Fe particle required a large 
bulk C supersaturation. Shibuta et al. [2l| also found a 
large C solubility in small Ni clusters, in contradiction 



with the ab initio calculations of Zhang et al. 14\ 
showing that, at K, C segregation at the surface of 
a M38 cluster is preferred. Studying the growth of 
nanotubes on small (M48 and Niso) clusters by empirical 
MD, Zhao et al. [22| found the same mechanism as in 
our previous calculations on a semi infinite system (slab 
geometry) [23| based on a tight binding model. These 
differences are partly due to the interatomic interaction 
model used but we show below that the initial state 
of the catalyst is an important parameter. Moreover, 
the MD technique itself is not able to take into account 
the carbon chemical potential ((J-c) gradient that is the 
thermodynamic driving force for the growth. 
In this Letter, we adopt a new approach, using a tight 
binding model for Ni-C interactions coupled to grand 
canonical Monte Carlo (GCMC) simulations, to gain 
a new insight on the of nucleation of C nanotubes. 
Experimentally, fic is imposed by the thermochemistry 
of the catalytic decomposition reaction of the carbon rich 
feedstock, which is the first step of all CVD processes. 
Introducing fj,c as a control variable of the computer 
simulations requires using a GCMC algorithm, with 
fixed volume, temperature, number of Ni atoms and 
Mc [HI- Such simulations are long and proceed by 
sequentially attempting random changes of the atomic 
configurations: they have to rely on a model in which 
the total energy is taken as a sum of local terms, to 
avoid recalculating it for the complete simulation box at 
each step. 

We developed and carefully tested [IE IIH such a 
model for Ni and C, in a tight binding framework. 
The total energy is a sum of local terms: an empirical 
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repulsive one and a band structure one including s and 
p electrons of C and d electrons of Ni. Local densities 
of electronic states (LDOS) are calculated using the 
recursion method. To keep the model as simple and fast 
to compute as possible, we neglect the Ni s electrons 
and calculate only the first four moments of the LDOS. 
The energy of each atom therefore depends only on the 
positions and chemical identities of its first and second 
neighbors, as defined by a cut off distance set at 3.20 A . 
In the GCMC algorithm used here, all atoms are allowed 
to move while only C atoms are tentatively added or 
deleted in a volume close (3 A above or below) to 
the surface of the Ni particle. Moreover, the lowest 
part (typically 40 %) of the particle is excluded for 
attempted insertions, in order to simulate a supported 
cluster and to avoid a possible encapsulation of the 
particle. More details on this algorithm that closely 
follows the thermodynamic conditions of CVD synthesis 
can be found in [26j. To analyze the structures, we 
distinguish three types of C atoms, according to the 
type (C or Ni) and length {dec or dcNi) of their nearest 
neighbor bonds. Outer C atoms, that ultimately form a 
nanotube cap, are defined as C atoms with at least one 
C neighbor at dec < 1-70 A . The remaining C atoms 
are subdivided into surface or bulk atoms, according to 
their number of Ni neighbors, defined by dcNi < 2.30 
A : those with less than 5 Ni neighbors are defined as 
surface C atoms, those with more than 5 neighbors as 
bulk C atoms. Note that these so called bulk atoms are 
essentially located in subsurface sites, but a diffusion 
towards interstitial sites buried deeper in the system is 
sometimes observed. 

The adsorption isotherms calculated at 1000 K, a 
typical temperature for CVD processes, are displayed in 
Fig. Q] Two clusters of 201 Ni atoms were considered 
: a crystalline one, in its equilibrium Wulff shape, and 
a disordered one, resulting from the quench of a liquid 
droplet. The fic values of interest for C adsorption 
are close to the energy of formation of graphene or 
diamond (-7.41 and -7.39 eV/atom in our model) and 
above. Adsorption begins at a lower fic value (-8.00 
eV/atom) on the FCC cluster than on the disordered 
one (-7.50 eV/atom) because the semi octahedral sites 
of the (100) facets of the former are very favorable. 
Above these values, and below a critical value \i* c that 
depends on the state of the Ni cluster, an increasingly 
large number of C atoms is adsorbed on the cluster 
surface or in subsurface interstitial sites. This n c value 
at which C atoms begin to form structures outside the 
cluster lies between -6.50 and -6.00 eV/atom for the 
crystalline cluster and between -6.00 and -5.75 eV/atom 
for the disordered one: below [L c the number of outer C 
atoms is zero, while a very steep increase is noticed at 
the transition. Around [i* c and beyond, the number of 
bulk and surface C atoms reaches a constant equilibrium 
value. We note that the number of bulk C atoms is 
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FIG. 1: Adsorption isotherms calculated on crystalline 
(black) or disordered (red) N1201 atoms clusters. Open circles 
and up pointing triangles correspond to C atoms dissolved in 
the bulk of the cluster; open squares and down pointing tri- 
angles to surface C atoms; full circles and open diamonds to 
outer C atoms. See text for the definition of these three kind 
of species. The arrows indicate that, beyond fic = -6.50 or 
-6.00 eV/at., the number of outer C atoms keeps increasing, 
while the other C species tend to equilibrium values. 



larger in the disordered Ni cluster (11±3) than in the 
crystalline one (6±2). This difference, observed on a 
nanosized object, is qualitatively the same as what is 
known for bulk phase diagrams: at constant tempera- 
ture, a liquid phase generally exhibits a larger solubility 
than the corresponding crystalline one. On the contrary, 
within the statistical error bars, the number of surface 
C atoms (34±3) does not depend on the state of the 
Ni particle. Referred to the 81 surface Ni atoms of the 
active part of the crystalline cluster, this corresponds 
to a critical surface C concentration of about 30 %. 
Canonical Monte Carlo calculations at 1000 K, starting 
with a fixed number of C atoms randomly deposited 
on a Ni(lll) surface, show that this critical surface 
concentration is large enough to allow some C atoms to 
interact with C neighbors and form small chains (with 3 
to 5 C atoms) that are stable over relatively long 'time' 
scales. Under grand canonical conditions, these chains 
act as nucleation sites onto which incoming C atoms 
attach. When the catalyst particle is disordered, the 
diffusion of C atoms towards subsurface sites is easier 
and the stability of these small chain embryos is lower, 
explaining the slightly larger fi* c value for disordered 
clusters. 

We now focus on the outer C atoms. Figure [5] presents 
a series of snapshots taken at different stages of the 
nucleation of a nanotube cap on a small disordered 
cluster of 50 Ni atoms, obtained by GCMC calculations 
at 1000 K and fic — -5.25 eV/atom. When the critical 
number of surface C atoms is reached, chains are formed, 
crawl on the surface and eventually cross each other. At 
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FIG. 2: GCMC calculations on a M50 cluster at 1000 K. Ni atoms are orange, outer C atoms are black, surface or bulk C 
atoms are green. Left panel: Final configuration obtained for fic = -6.50 eV : no growth is observed. Central panel : different 
stages of the nucleation of a C cap at fj-c = -5.25 eV/at.. (a) atoms or dimers adsorbed; (b) chains forming; (c) chains crossing; 
(d) cap formed; (e) cap lifts off. Right panel : final configuration obtained for nc = -4.50 eV resulting in the formation of an 
amorphous C layer. 



their intersections, three fold coordinated C atoms act as 
nucleation centers for C sp 2 structures that develop on 
the surface. These sp 2 C atoms interact weakly with the 
underlying Ni atoms and can detach from the surface. 
The energy of adhesion of a perfect graphene layer in 
epitaxy on a Ni( 111) surface has indeed been shown 
to be close to zero [H, [23| . This growth pattern was 
already observed on flat Ni slabs [23j. The difference 
is that the curvature necessary to form a nanotube 
cap is now provided by the curvature of the small Ni 
cluster, meaning that the diameter of the nanotube cap 
corresponds to the local curvature of the Ni particle, 
at the moment of the nucleation, as observed in the in 
situ TEM observations [lo| . In addition, the systematic 
analysis of the role of nc shows that there is an optimal 
He range for nucleating nanotubes. The existence of the 
lower threshold is in qualitative agreement with 

the conclusions of Lolli et al. [5], who showed that a 
lowering of the C surface fugacity hinders the nucleation 
of C species on the metal surface. Beyond the upper 
threshold (roughly, \ic > —5.00 eV/atom), the fast 
growth of an outer C structure leads to the formation of 
amorphous carbon (ultimately carbon fibers), similar to 
those observed by Helveg et al. 

We now understand the role of the catalyst, beyond 
its primary, truly catalytic, activity in the gas feed 
decomposition reaction. The above analysis of the 
partition of adsorbed C atoms into bulk, surface or 
outer sites shows that the first role of the metal catalyst 
is to confine atomic C on, or close to the surface until 
a critical concentration is reached. The local energies 
distributions plotted in Fig. [3] tell more. The energy 
distribution of Ni atoms is centered around -4.0 eV and 
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FIG. 3: Energy distribution of each species, calculated on 
Ni2oi crystalline (lower panel) or disordered (upper panel). 
Note the bimodal shape of the outer C atoms distribution. 



that of outer C atoms is bimodal: the peak around -5.8 
eV corresponds to C atoms with 2 C neighbors, while 
the peak around -7.2 eV corresponds to the more stable 
threefold coordinated C atoms. Using in situ XPS 
spectroscopy Hofmann et al. 10] evidenced the presence 
of chemisorbed C close to the surface of the Fe catalyst, 
stable during ~ 150 s, prior to the formation of graphitic 
C. We also observe this and note that the energy of 
individual C atoms adsorbed on the surface matches 
that of the twofold coordinated outer C atoms: this is 
the second role of the metal catalyst. As illustrated 



4 



in Fig. [2 it makes the transition from isolated atoms 
or dimers adsorbed on the surface to chains of sp C 
possible with a low energy cost, as soon as the critical 
concentration is reached. This transition is a continuous 
one, without activation energy needed. As also shown 
in [15J, |25| subsurface C atoms are only slightly more 
stable than surface ones, otherwise the tendency to form 
a carbide would prevent the surface segregation of C. 
Deck and Vecchio have shown experimentally that 
successful catalysts have C solubility limits between 0.5 
wt.% and 1.5 wt.% C. The energy pattern observed is 
compatible with such a limited solubility range and is 
the key for nucleating C nanotubes. 
As shown in Fig. [4j the bulk or subsurface carbon 
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FIG. 4: Final state of disordered (a) or initially crystalline 
(b) M201 catalyst particles, after formation of an outer C 
layer that has been removed for legibility. Note that the lat- 
ter remains essentially crystalline, except for the outermost 
layer. Green atoms are surface C atoms, blues ones are bulk 
C atoms. See text for definition of these two kinds of species. 

concentration depends on the state (crystalline or 
disordered) of the catalyst. The initially crystalline 
cluster retains a large degree of crystallinity, but the 
observed preference for incorporating C atoms in in- 
terstitial sites located below edge or surface Ni atoms 
induces a strong local disorder. As shown in (23[, a large 
concentration of C in the outermost Ni layers induces a 
surface melting. Depending on the initial state of the 
particle, the temperature and the chemical potential, a 
molten surface layer may exist, making a fast diffusion 
of both Ni and C possible and possibly explaining the 
reshaping of the catalyst particle observed by in situ 
TEM experiments 0, [l(| ■ We thus reconcile conflicting 
experimental observations : some TEM observations 
(Zhu et at [g], for example) show clearly crystalline 
particles, while others (e.g.: Harutyunyan et al. |28j |) 
and time resolved reflectivity measurements by Puretzky 
et al. support the presence of a molten or partially 
molten catalytic particle. 

In this Letter we have shown that an optimal \xc window 
exists to nucleate nanotube caps whose curvature match 
the local curvature of the catalyst particle. The nucle- 
ation is triggered when the concentration of C atoms 
adsorbed on the surface is large enough, about 30 % in 



our calculations. The role of the catalyst is to confine 
C atoms on or close to the surface, and to make them 
reach this critical concentration. Moreover, the energies 
of isolated C atoms on the surface and of C chains 
adsorbed are in the same range, making a continuous 
transition possible. The C concentration in bulk and 
subsurface layers depends on the state (crystalline or 
disordered) of the catalyst. Beyond their agreement with 
most experimental observations, these results emphasize 
the importance of the limited C solubility in Ni and 
give an insight at the atomic level on the structure of 
the tube/catalyst interface. They represent an essential 
step towards a better control of the structure (diameter, 
chirality) of the tubes formed and should help develop 
better catalysts. 
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